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Abstract 

The  direct  numerical  simulation  (DNS)  data  of  a 
turbulent,  reacting  boundary  layer  is  used  to  study  the 
turbulence-chemistry  interaction  and  the  scaling  of  the 
temperature  fluctuations.  We  find  that  there  is  a  feed¬ 
back  between  the  turbulence  and  the  chemical  reac¬ 
tions.  Temperature  fluctuations  increase  the  reaction 
rates.  Endothermic  reactions  reduce  the  magnitude  of 
the  temperature  fluctuations.  However,  the  turbulent 
boundary  layer  is  self-sustained  and  the  magnitude  of 
the  temperature  fluctuations  is  maintained.  Exother¬ 
mic  reactions  increase  the  magnitude  of  the  tempera¬ 
ture  fluctuations.  The  reactions  cause  localized  com¬ 
pressions  and  expansions  that  feed  the  turbulent  ki¬ 
netic  energy.  The  evolution  equations  for  the  turbu¬ 
lent  energy  show  that  the  feedback  takes  place  through 
the  transport  and  production  terms.  The  data  reveals 
a  physically  consistent  relation  between  the  tempera¬ 
ture  variance  and  the  mean  flow  parameters.  These 
results  are  consistent  with  those  previously  found  in 
reacting  isotropic  turbulence.1,2 

Introduction 

The  boundary  layers  that  are  formed  on  hypersonic 
vehicles  are  hot,  chemically  reacting,  and  turbulent. 
Currently,  the  boundary  layer  on  realistic  hypersonic 
vehicles  is  simulated  either  assuming  that  the  bound¬ 
ary  layer  is  laminar  or  using  simple  models  that  have 
not  been  calibrated  for  hypersonic  applications.  Gen¬ 
erally,  the  calibration  of  turbulent  models  has  been 
done  using  DNS  databases  of  incompressible  flows  or 
using  perfect  gas  wind  tunnel  data.  If  we  were  able 
to  perform  more  accurate  simulations  of  hypersonic 
flows,  we  may  find  a  different  chemical  composition  of 
the  gas,  and  different  heating-rates  than  those  that  are 
currently  predicted. 

In  the  flows  of  interest  the  magnitude  of  the  temper¬ 
ature  fluctuations  is  very  large  due  to  the  high  energy 
content  that  is  present.  Furthermore,  the  chemical 
reaction  is  a  highly  non-linear  function  of  the  temper¬ 
ature.  Therefore,  temperature  fluctuations  may  result 
in  large  variations  in  the  reaction  rates.  In  this  paper 
we  study  the  turbulence-chemistry  interaction  in  a  tur¬ 
bulent  boundary  layer  and  seek  a  scaling  of  the  tem¬ 
perature  fluctuations  based  on  the  resolved  mean  flow 


variables.  The  aim  is  to  devise  a  model  for  the  tem¬ 
perature  fluctuations  that  can  be  used  in  the  context 
of  Reynolds-averaged  Navier-Stokes  or  large-eddy  sim¬ 
ulations  to  predict  accurate  heating-rates  and  product 
formation. 

Following  our  previous  work,3  we  use  the  DNS 
database  of  a  boundary  layer  at  Reg  =  7000,  edge- 
conditions  of  Me  =  4.0,  pe  =  0.5  kg/m3,  Te  =  5000 
K,  and  wall-temperature  Tw  =  5000  K.  These  condi¬ 
tions  represent  the  boundary  layer  on  a  26°  wedge  at  a 
Mach  number  of  20  and  20  km  altitude.  We  use  a  sin¬ 
gle  model  reaction,  SI  +  M  #  S2  +  M,  where  species 
SI  and  S2  have  the  same  molecular  weight  and  number 
of  degrees  of  freedom.  In  this  way,  the  gas  constant  is 
not  a  function  of  the  chemical  composition  of  the  gas 
and  changes  in  pressure  are  only  caused  by  density 
and  temperature  variations.  The  reaction  rates  corre¬ 
spond  to  oxygen  dissociation,  and  SI  and  S2  represent 
molecules  and  atoms,  respectively.  Thus,  production 
of  species  SI  indicates  an  exothermic  reaction.  The 
reaction  rate  and  equilibrium  constant  expressions  are 
taken  from  Gupta  et  al.4 

The  paper  is  organized  as  follows.  We  first  present 
a  review  of  our  previous  work  by  introducing  the  non- 
dimensional  parameters  that  govern  the  turbulence- 
chemistry  interaction  and  the  scaling  functions  that 
were  found  using  DNS  data  of  isotropic  turbulence  at 
conditions  typical  of  hypersonic  boundary  layers.  In  a 
turbulent  boundary  layer,  the  energy  transfer  mecha¬ 
nisms  are  more  complex  than  those  found  in  isotropic 
turbulence.  Thus,  we  present  the  temporal  evolution 
equations  for  the  turbulent  kinetic  energy  and  the  vari¬ 
ance  of  the  temperature.  These  equations  are  used  in 
the  analysis  of  the  turbulence-chemistry  interaction. 
The  results  are  divided  in  two  sections.  First,  we 
present  the  characteristics  of  the  mean  flow  and  the 
effect  of  turbulent  fluctuations  on  the  mean  flow  mo¬ 
tion.  We  then  discuss  the  turbulent  energy  transfer 
and  the  scaling  of  the  temperature  fluctuations.  Con¬ 
clusions  are  given  in  the  last  section. 
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Background 

The  nondimensional  parameters  governing  the  tur¬ 
bulence  chemistry  interaction  are1  the  turbulent  Mach 
number,  the  Reynolds  number,  the  Damkohler  num¬ 
ber,  and  the  relative  heat  release,  namely 
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where  q  =  (it'u^)1/2  is  the  rms  magnitude  of  the  fluc¬ 
tuation  velocity;  a  is  the  speed  of  sound;  t*  =  k/e 
is  the  turbulent  time  scale,  where  k  and  e  are  the 
turbulent  kinetic  energy  and  dissipation,  respectively; 
Tc  =  pl  ((w|1)1/2Ke,)  is  the  chemical  time  scale  where 
w  is  the  source  term,  and  Ke<?  is  the  equilibrium  con¬ 
stant;  A  =  qrt  is  the  characteristic  turbulent  length 
scale;  and  A h°  is  the  heat  of  the  reaction. 

In  our  previous  work,1,2  we  used  DNS  to  perform  a 
fundamental  study  of  isotropic  turbulence  interacting 
with  finite-rate  chemical  reactions  at  conditions  typ¬ 
ical  of  a  hypersonic  boundary  layer.  We  found  that 
the  turbulent  motion  is  fed  from  the  energy  provided 
by  the  exothermic  reactions,  while  the  reaction  rate  is 
increased  by  the  turbulent  temperature  fluctuations. 
This  is  a  feedback  process  that  takes  place  through  the 
pressure-strain  term  in  the  Reynolds  stress  equation. 
The  feedback  is  negative  for  endothermic  reactions, 
resulting  in  a  reduction  in  the  turbulent  motion.  The 
DNS  database  showed  that  the  temperature  fluctua¬ 
tions  can  be  expressed  as  a  function  of  the  governing 
parameters.  For  endothermic  reactions,  the  tempera¬ 
ture  fluctuations  scale  linearly  with  M2.  Whereas  for 
exothermic  reactions,  the  temperature  fluctuations  are 
enhanced  and  can  be  expressed  as 


T'Rms/(T )  =  A{£xh°  \/lE)B 
=  A(Ah°  MtDa)B 


(2) 


where  A  and  B  are  constants  that  depend  on  the  spe¬ 
cific  reaction  and  A,  also  known  as  as  the  Taylor  mi¬ 
croscale,  represents  the  distance  traveled  by  a  fluid 
particle  moving  at  the  speed  of  the  turbulent  inten¬ 
sity.  The  expansion  length  Ie  is  defined  as  arc,  which 
is  the  distance  traveled  by  acoustic  radiation  from  the 
chemistry-induced  temperature  fluctuations.  There¬ 
fore,  A /Ie  represents  the  ratio  of  the  characteristic 
distance  traveled  by  a  typical  particle  of  fluid  to  the 
characteristic  distance  traveled  by  the  acoustic  radia¬ 
tion. 

Another  way  to  understand  this  ratio  of  length 
scales  is  to  consider  the  variation  of  the  strength  of  the 


chemistry-turbulence  interaction.  A  positive  tempera¬ 
ture  fluctuation  increases  the  reaction  rate,  making  the 
reaction  occur  more  quickly,  which  releases  more  heat, 
further  increasing  the  temperature.  However,  the  feed¬ 
back  process  can  be  weakened  by  delocalization  of  the 
interaction  through  turbulent  motion  and  motion  gen¬ 
erated  by  the  local  pressure  fluctuations  (also  caused 
by  the  interaction).  Therefore,  in  a  simplified  case 
where  the  reaction  rate  is  held  constant,  the  strength 
of  the  interaction  varies  like  (u’)  -^-r-  =  —  =  In 
general,  the  interaction  strength  also  varies  with  the 
reaction  rate,  or  in  non-dimensional  terms,  with  the 
Damkohler  number.  Thus  with  this  argument,  we  ob¬ 
tain  the  result  shown  in  eq.  (2). 

Under  the  conditions  chosen  for  our  calculations  of 
isotropic  turbulence,  A /Ie  was  always  less  than  one. 
As  A /Ie  approaches  one,  FrMS/(T)  becomes  large, 
indicating  a  strong  turbulence-chemistry  interaction. 
This  occurs  when  the  fluid  travels  a  similar  distance 
as  the  acoustic  radiation  induced  by  the  temperature 
fluctuations.  If  A /Ie  were  larger  than  one,  the  interac¬ 
tion  would  be  expected  to  weaken  because  the  turbu¬ 
lent  motion  would  outrun  the  acoustic  waves  produced 
by  the  interaction,  and  the  feedback  process  would  be 
diminished.  Also,  as  A /Ie  approaches  zero  the  pressure 
waves  outrun  the  fluid  motion  and  the  interaction  is 
weak.  Thus,  the  interaction  weakens  when  A /Ie  de¬ 
parts  from  unity.  In  addition,  TrMS/(T)  is  affected 
by  the  heat  released  to  the  flow,1  and  the  length  ra¬ 
tio  must  be  modulated  by  A h°  to  give  an  appropri¬ 
ate  relation  for  the  standard  deviation.  When  A/Ie 
is  greater  than  one,  we  would  not  expect  this  fit  to 
be  valid  because  it  predicts  a  further  strengthening  of 
the  interaction.  The  Reynolds  number  did  not  have 
a  significant  effect  for  the  range  of  conditions  that  we 
considered.  A  more  detailed  discussion  about  this  scal¬ 
ing  for  the  temperature  fluctuations  can  be  found  in 
Martin  &  Candler.2 


Turbulence  mechanisms 

The  chemical  reactions  act  as  energy  sources  within 
the  turbulent  boundary  layer.  Thus,  we  must  address 
the  energy  exchange  between  the  turbulence  and  the 
chemical  reactions.  There  are  four  energy  exchange 
mechanisms  that  take  place  in  turbulent  boundary  lay¬ 
ers:  transport,  production,  dissipation  and  diffusion  of 
turbulence.  The  budget  equation  for  the  turbulent  ki¬ 
netic  energy  is 

d  d 

—(pk)  +  w—(p  k )  =  Pk  +  Tk  +  Bt  +Ud  +  (j>dif  +  (f>dis 
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and  Pfc  is  the  production  due  to  the  mean  gradients, 
Tfc  is  the  redistribution  of  turbulent  kinetic  energy,  IR 
is  the  pressure  diffusion,  and  11^  is  the  pressure  dila¬ 
tion,  (j>dif  is  the  viscous  diffusion,  (f>diS  is  the  viscous 
dissipation.  There  are  other  terms  that  appear  in  the 
equation  due  to  the  Favre  averaging,  however  these 
terms  are  negligible.  Note  that  v!  and  u"  represent 
fluctuations  with  respect  to  the  Reynolds  and  Favre 
averages  of  u,  respectively. 

To  study  the  turbulent  internal  energy  we  use  the 
evolution  equation  for  the  temperature  variance,  which 
is  given  by 


z  /  5 


Figure  1.  Average  (a)  temperature  and  (b)  mass 
fraction  of  species  SI. 
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and  Pt  is  the  production  of  temperature  fluctuations 
due  to  the  chemical  reactions  and  h°  is  the  heat  of 
formation  of  species  s;  Tt  is  the  transport  of  temper¬ 
ature  fluctuations  due  to  the  temperature-dilatation 
correlation;  ipvdiS  represents  the  viscous  dissipation  of 
kinetic  energy  into  internal  energy;  ptdif  is  the  dif¬ 
fusion  of  temperature  fluctuations  due  to  heat  con¬ 
duction  qj  =  —ndT/dxj]  and  ipsdif  is  the  redistribu¬ 
tion  of  temperature  fluctuations  due  to  the  diffusion 
of  species,  where  D  is  the  species  diffusivity  and  cs  is 
the  mass  fraction  of  species  s.  The  last  term,  ipsdif  is 


negligible  since  for  our  model  reaction  the  gradients  of 
species  SI  and  S2  are  equal  in  magnitude  and  opposite 
in  sign. 


Flow  characteristics 

The  species  mass  fractions  are  initialized  to  their 
equilibrium  values  at  the  averaged  temperature.  This 
initial  state  is  not  physical.  Thus,  when  we  turn  on 
the  reactions  the  flow  undergoes  a  transient.  Figure  la 
shows  the  average  temperature  for  the  initial  condition 
and  after  the  transient.  The  gas  near  the  wall  is  cooled 
considerably.  The  temporal  evolution  equation  for  the 
average  temperature  shows  that  the  mechanisms  caus¬ 
ing  this  effect  are  the  chemical  source  term  and  the 
thermal  diffusion.  Figure  lb,  however,  shows  that  the 
mass  fraction  of  species  SI  increases  through  the  tran¬ 
sient.  Production  of  SI  indicates  exothermic  reactions, 
which  would  result  in  an  increase  of  the  average  tem¬ 
perature.  This  result  can  be  explain  by  looking  at  the 
budget  of  the  terms  in  the  evolution  equation  for  the 
density  of  species  SI,  which  shows  that,  locally,  species 
SI  diffuses  at  a  faster  rate  than  the  source  term  is  de¬ 
stroying  it.  After  the  initial  transient,  the  mean  flow 
is  in  a  quasi-equilibrium  state.  Although  the  chemi¬ 
cal  composition  in  typical  hypersonic  boundary  layers 
is  not  in  chemical  equilibrium,  these  flow  conditions 
serve  to  isolate  the  effect  of  turbulent  fluctuations  in 
the  turbulence-chemistry  interaction. 
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Figure  2.  (a)  Average  production  of  species  SI  nor¬ 
malized  by  pwuT/zT ;  (b)  budget  for  the  time  evolution 
of  the  average  temperature  normalized  by  u\jcv  zT. 


Figure  3.  Normalized  magnitude  of  the  fluctuations 
in  the  production  of  species  SI. 


Figure  2a  shows  the  average  production  of  SI. 
There  is  nearly  no  net  chemical  production  away  from 
the  wall.  At  the  wall  ws  is  negative,  indicating  that  de¬ 
struction  of  SI  or  endothermic  reactions  are  dominant. 
As  mentioned  above,  the  mean  flow  is  in  chemical  equi¬ 
librium,  thus  the  average  temperature  should  not  de¬ 
crease  due  to  the  endothermic  reactions.  Figure  2b 
shows  the  budget  for  the  terms  in  the  evolution  equa¬ 
tion  for  ( T)t  .  This  equation  is  similar  to  the  equation 
for  (T'T')t  but  without  the  factors  of  2 T'  in  the  terms 
on  the  right-hand  side  of  eq.  (4).  We  observe  that 
the  effect  of  the  endothermic  reactions  is  to  decrease 
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Figure  4.  Normalized  magnitude  of  the  fluctuations 
in  (a)  temperature  and  (b)  mass  fraction  of  species  SI. 

the  temperature.  Also  the  effect  of  the  thermal  diffu¬ 
sion  is  to  cool  the  gas  near  the  wall.  However,  these 
mechanisms  are  balanced  by  the  large  dissipation  of 
kinetic  energy  into  heat,  which  is  represented  by  the 
viscous  dissipation. 

We  now  consider  the  effect  of  the  turbulent  fluc¬ 
tuations  in  the  mean  flow.  Figure  3  shows  the  rms 
of  the  fluctuating  source  term  normalized  by  its  abso¬ 
lute  magnitude.  The  magnitude  of  the  fluctuations  are 
higher  than  60%  of  the  average  production  at  the  wall 
and  100%  elsewhere,  indicating  that  the  production 
of  species  is  mainly  due  to  fluctuations  in  the  source 
term.  The  temperature  fluctuations  drive  the  chemical 
source  term  and  since  the  mean  flow  is  in  equilibrium, 
no  effect  is  observed  in  the  average  chemical  compo¬ 
sition  of  the  gas,  namely  (csi)  equals  its  equilibrium 
value  based  on  the  averaged  temperature.  However,  we 
should  be  able  to  observe  an  effect  on  the  magnitude 
of  the  fluctuating  csi-  Figure  4a  shows  the  magnitude 
of  the  temperature  fluctuations  after  the  initial  tran¬ 
sient.  The  maximum  is  roughly  5%  of  the  average  and 
takes  place  in  the  viscous  sublayer.  A  second  peak  of 
nearly  the  same  magnitude  develops  in  the  logarithmic 
region.  These  locations  are  marked  with  symbols.  For 
a  given  temperature  there  is  a  corresponding  equilib¬ 
rium  composition  of  the  gas.4  Table  1  gives  the  values 
of  csi  evaluated  at  (T)  +  T'RMS  and  (T)  —  T'RMS  and 
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□ 

X 

Csi  «T)  +T’RMS)  /{csi) 

25% 

24% 

csi  ((T)  —  T'rms )  /(cSi) 

31% 

31% 

Table  1 .  Equilibrium  mass  fraction  of  species  SI  eval¬ 
uated  at  the  maximum  and  minimum  local  tempera¬ 
tures  in  the  boundary  layer  and  relative  to  its  mean. 


normalized  by  (csi)  (which  equals  csi  evaluated  at 
(T)).  We  see  that  5%  fluctuation  in  the  mean  tem¬ 
perature  can  cause  a  31%  difference  in  the  chemical 
composition.  Clearly,  the  equilibrium  composition  is 
a  very  strong  function  of  the  temperature.  This  is  re¬ 
flected  in  the  magnitude  of  c'si  within  the  boundary 
layer,  which  is  shown  in  Fig.  4b.  We  see  that  varia¬ 
tions  in  the  equilibrium  mass  fraction  are  nearly  that 
produced  by  the  temperature  fluctuations. 


Turbulence-chemistry  interaction 

In  this  section  we  study  the  turbulence-chemistry 
interaction  and  the  temperature  fluctuation  scaling  by 
studying  the  budget  of  the  terms  in  the  evolution  equa¬ 
tions  for  the  temperature  variance  and  the  turbulent 
kinetic  energy.  Figure  5  shows  the  effect  of  the  source 
term  on  the  temperature  fluctuation  variance.  The 
symbols  are  shown  for  future  reference.  The  source 
term  has  been  decomposed  into  the  forward  and  back¬ 
ward  components, 


w  si  = 


vr  i  PSifPSi  .  PS2  \ 

S1  f  Msi  VMsi  MS2) 

,  ,  PS2  /  PSl  .  PS2  \ 

+  Msi  h  Ms^MsT  +  Ms^J 

=  Wf  +  Wb, 


(6) 


where  ps  and  Ms  are  the  density  and  molecular  weight 
of  species  s,  kf  and  kb  are  the  forward  and  back¬ 
ward  reactions  rates,  respectively,  and  Wf  and  Wb  rep¬ 
resent  the  destruction  and  production  of  species  SI. 
Thus,  Wf  reduces  the  magnitude  of  the  temperature 
fluctuations  as  a  result  of  endothermic  reactions  and 
Wb  has  the  opposite  effect.  Figure  5  also  shows  that 
both  exothermic  and  endothermic  reactions  take  place 
within  the  boundary  layer.  Figure  6a  shows  the  budget 
of  ( T'T ')  in  the  near  wall  region.  At  the  wall  the  chem¬ 
ical  production  and  thermal  diffusion  are  balanced 


Figure  5.  Contribution  of  the  chemical  production 
to  ( T'T')t  ■  The  terms  are  normalized  using  u?Tjcv  zT. 


Figure  6.  Budget  of  the  terms  in  the  evolution 
equation  for  (T'T')t  (a)  near  the  wall  and  (b)  across 
the  boundary  layer.  The  terms  are  normalized  using 
u^/cv  zT. 

by  the  viscous  dissipation  and  transport  mechanisms. 
The  species  diffusivity  is  negligible  since  the  gradients 
of  species  mass  fraction  are  equal  in  magnitude  and 
opposite  in  sign.  Figure  6b  shows  the  same  budget 
throughout  the  entire  boundary  layer.  The  chemi¬ 
cal  production  and  transport  mechanisms  enhance  the 
temperature  fluctuations,  however  these  terms  are  bal¬ 
anced  by  the  thermal  diffusion.  The  species  diffusion 
term  is  negligible  and  not  shown. 

Figure  7a  shows  the  budget  of  the  turbulent  kinetic 
energy  in  the  near  wall  region.  The  turbulent  kinetic 
energy  is  nearly  in  equilibrium.  Figure  7b  shows 
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Figure  7.  Budget  of  the  terms  in  the  evolution 
equation  for  (pk)t  (a)  near  the  wall  and  (b)  across 
the  boundary  layer.  The  terms  are  normalized  using 
pwu\/zT. 


Figure  8.  (a)  Temporal  evolution  of  the  wall- friction 
velocity  and  (b)  turbulent  kinetic  energy  profile  nor¬ 
malized  using  pwv%- 


the  same  budget  across  the  entire  boundary  layer.  The 
production  term  is  dominant  causing  a  net  produc¬ 
tion  of  turbulent  kinetic  energy.  The  production  and 
transport  terms  have  the  same  oscillatory  behavior  in 
0.1  <  z/S  <  0.3.  This  behavior  is  similar  to  that  of  the 
transport  term  in  the  equation  for  (T'T')t,  as  shown 
in  Fig.  6b.  Figure  8a  plots  the  temporal  evolution 
of  the  friction  velocity,  uT.  After  the  initial  transient, 
the  oscillations  in  uT  indicate  the  presence  of  a  phys¬ 
ical  mechanism  that  is  damping  and  replenishing  uT. 
Figure  8b  plots  the  turbulent  kinetic  energy  across  the 
boundary  layer.  Again  the  oscillatory  behavior  is  ob¬ 
served  in  0.1  <  z/S  <  0.3,  as  in  the  transport  and 
production  mechanisms  for  the  evolution  of  (T'T')t 
and  (TKE)t.  As  in  isotropic  turbulence,1  this  result 
suggests  that  a  feedback  mechanism  between  the  tur¬ 
bulent  motion  and  the  chemical  reactions  takes  place. 

The  effect  of  endothermic  reactions  is  to  damp  the 
temperature  fluctuations.  However  for  the  conditions 
chosen,  the  turbulent  kinetic  energy  in  the  boundary 
layer  is  self-sustained.  Therefore,  the  turbulent  tem¬ 
perature  fluctuations  are  maintained.  The  exother¬ 
mic  reactions,  however,  act  as  heat  sources  increas¬ 
ing  the  temperature  fluctuations.  Both  endother¬ 
mic  and  reactions  cause  localized  compressions  and 


exothermic  expansions  within  the  boundary  layer.  The 
turbulent  kinetic  energy  feeds  from  the  thermal  en¬ 
ergy.  The  results  suggest  that  this  feedback  mecha¬ 
nism  takes  place  through  the  transport  term  in  the 
equation  for  ( T'T')t  and  through  the  production  and 
transport  mechanisms  in  the  equation  for  (TKE)t. 
The  transport  term  for  (T'T')t  includes  the  correla¬ 
tion  between  temperature  fluctuations  and  dilatation, 
which  explains  the  energy  transfer  between  the  turbu¬ 
lent  kinetic  energy  and  the  internal  energy  through  the 
chemically  induced  expansions  and  compressions. 

Figure  9  plots  the  magnitude  of  the  temperature 
fluctuations  versus  the  heat  ratio,  eq.  (2),  in  the  vis¬ 
cous  sublayer  and  near  z/S  =  0.5.  (  These  boundary 
layer  locations  are  marked  with  symbols  in  Fig.  5.) 
The  data  includes  several  time  steps  during  the  sim¬ 
ulation.  The  power  law  relation  between  temperature 
fluctuations  and  the  heat  ratio  resembles  that  found 
in  isotropic  turbulence.  The  data  collapse  very  well 
in  the  viscous  sublayer,  where  the  turbulence  scales 
are  nearly  isotropic.  Near  the  boundary  layer  edge 
however,  the  agreement  is  not  so  good.  This  could  be 
because  near  the  boundary  layer  the  numerical  reso¬ 
lution  is  not  as  good,  since  the  computational  grid  is 
exponentially  stretched.  This  can  be  seen  in  the  lower 
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Figure  9.  Temperature  fluctuation  versus  heat  ratio. 


population  of  symbols  for  the  data  near  the  bound¬ 
ary  layer.  It  is  also  possible  that  we  have  not  gather 
sufficient  statistics  in  time  so  as  to  get  a  large  enough 
sample  of  data.  A  Reynolds  number  dependence  is  also 
possible,  since  the  peak  in  the  logarithmic  region  of  the 
turbulent  profiles  is  highly  dependent  on  the  Reynolds 
number.  None  of  these  reasons  or  others  can  be  proven 
true  with  the  small  dataset  that  we  are  working  with. 

Conclusions 

In  this  paper,  we  have  analyzed  the  direct  numeri¬ 
cal  simulation  database  of  a  chemically  reacting,  tur¬ 
bulent  boundary  layer  at  Mach  4,  Reg  =  7000,  and 
isothermal  wall-temperature  condition.  The  budget  of 
the  terms  in  the  evolution  equation  for  the  tempera¬ 
ture  show  that  the  mean  flow  is  in  chemical  equilib¬ 
rium.  However,  the  temperature  fluctuations  increase 
the  reaction  rates  and  therefore  the  chemistry  is  ac¬ 
tive  throughout  the  boundary  layer.  The  magnitude  of 
the  fluctuations  in  temperature  and  species  mass  frac¬ 
tion  show  that  a  5%  fluctuation  in  temperature  can 
cause  a  30%  fluctuation  in  the  species  mass  fraction. 
If  the  boundary  layer  is  in  chemical  non-equilibrium, 
such  large  fluctuations  in  the  species  mass  fraction  may 
have  a  significant  effect  in  the  average  composition  of 
the  gas  and  therefore  in  the  heating  rates  to  the  wall. 
Additional  simulations  are  required  to  verify  this. 

The  evolution  equation  for  the  temperature  fluctua¬ 
tions  show  that  endothermic  reactions  reduce  the  tem¬ 
perature  fluctuations,  whereas  exothermic  reactions 
act  as  heat  sources  that  increase  T'RMS.  The  budget 
of  turbulent  kinetic  energy  shows  that  the  production 
and  transport  mechanisms  are  enhanced  by  the  pres¬ 
ence  of  reactions.  The  oscillations  found  in  the  evo¬ 
lution  of  the  wall-friction  velocity  and  in  the  profile 
of  turbulent  kinetic  energy  indicate  the  presence  of  a 
physical  mechanism  that  is  damping  and  replenishing 
the  turbulent  kinetic  energy. 

The  effect  of  endothermic  reactions  is  to  damp  the 
temperature  fluctuations.  However,  the  turbulent  ki¬ 


netic  energy  in  the  boundary  layer  is  self-sustained. 
Therefore,  the  turbulent  temperature  fluctuations  are 
maintained.  As  mentioned  above,  exothermic  reac¬ 
tions  increase  the  magnitude  of  the  temperature  fluc¬ 
tuations  and  the  temperature  fluctuations  increase  the 
reaction  rates.  Both  endothermic  and  exothermic  re¬ 
actions  cause  localized  compressions  and  expansions 
within  the  boundary  layer  that  feed  the  turbulent  ki¬ 
netic  energy.  The  results  suggest  that  this  feedback 
mechanism  takes  place  through  the  transport  term  in 
the  equation  for  {T'T')t  and  through  the  production 
and  transport  mechanisms  in  the  equation  for  ( pk)t ■ 
The  transport  term  for  (T'T')t  includes  the  correla¬ 
tion  between  temperature  fluctuations  and  dilatation, 
which  explains  the  energy  transfer  between  the  tur¬ 
bulent  kinetic  energy  and  the  internal  energy  through 
the  chemically  induced  expansions  and  compressions. 
This  result  is  consistent  with  that  found  using  the  DNS 
data  of  reacting  isotropic  turbulence1,2  and  with  the 
theoretical  work  of  Eschenroeder.5 

Also  as  in  isotropic  turbulence,  we  find  that  the 
temperature  variance  can  be  expressed  as  A  (AIi°\/Ie) 
Where  X/Ie  is  obtained  from  the  nondimensional  gov¬ 
erning  parameters  and  represents  the  ratio  of  the  char¬ 
acteristic  distance  traveled  by  a  fluid  particle  to  the 
characteristic  distance  of  acoustic  radiation.  Since 
Ah°X/lE  is  a  mean  flow  quantity,  it  could  potentially 
be  used  to  predict  the  temperature  fluctuations  in 
Reynolds-averaged  Navier-Stokes  or  large-eddy  simu¬ 
lations.  However,  further  testing  of  the  scaling  must 
be  performed  over  a  wider  range  of  conditions. 
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